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THE  EMPTYING  OF  A  GUN  TUBE  (U) 

ABSTRACT 

An  elementary  approximate  description  is  proposed  for  the  flow 
of  propellant  gases  from  a  gun  tube.  After  a  recapitulation  of  the 
discussion  of  the  Pidduck-Kent  flow,  there  is  presented  a  simple  flow 
which  is  conjectured  to  he,  with  the  unknown  correct  value  of  a 
certain  parameter,  an  accurate  asymptotic  representation  of  the  (.suit¬ 
ably  idealized)  succeeding  flow.  This  Binrple  flow  can  be  generalized 
to  the  emptying  of  an  arbitrary  vessel  through  an  orifice  and  is  in 
turn  a  generalization  of  the  quasi  stationary  treatment  of  the  emptying 
through  a  small  orifice.  In  the  gun  it  appears  to  be  a  good  approxi¬ 
mation  except  for  a  brief  interval  after  the  disturbance  of  the  flow 
by  the  release  of  pressure  at  the  muzzle.  An  example  is  given  in  which 
numerical  comparison  can  be  made  with  experiment  and  with  other  calcu¬ 
lations  . 
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The  flow  of  propellant  gases  from  a  gun  tube  becomes  Important 
in  a  number  of  problems ,  among  which  are  l)  the  estimation  of  recoil, 

2)  the  design  for  very  rapid  fire  in  which  the  breech  may  be  opened 
before  the  pressure  has  completely  decayed,  3)  the  design  of  a  muzzle 
brake  or  bore  evacuator.  We  are  presenting  here  a  simple  analytic 
treatment  of  the  emptying  of  the  tube  which  may  be  useful  for  (l)  and 
(2)  but  not,- in  many  cases,  for  such  purposes  as  (3)  since  It  provides 
a  good  approximation  to  the  flow  except  just  after  the  disturbance  of 
the  flow  hy  the  release  of  pressure  at  the  muzzle.  The  aim  Is  to  ob¬ 
tain  a  simple  realistic  treatment  which  can  be  used  with  confidence  In 
extreme  cases  and  over  long  intervals  of  time.  Existing  methods  appear 
to  he  sufficiently  accurate  for  most  needs  in  conventional  cases.  The 
flow  of  the  gas  from  the  tube  is  pictured  as  consisting  of  two  stages. 
The  first  stage  is  governed  completely  by  the  initial  condition,  the 
flow  attained  by  the  expanding  column  of  gas  at  the  Instant  the  base 
of  the  projectile  emerges  from  the  muzzle.  We  have  supposed  this 
portion  of  the  flow  to  he  the  flow  described  by  Love  and  Pidduck  [8]  and 
Kent  (5)  .  In  the  first  section  of  this  report  we  have  recapitulated, 
perhaps  unnecessarily,  the  discussion  of  that  flow  and  have,  we  hope, 
simplified  it  somewhat. 

The  second  stage,  which  we  have  called  the  expansion  wave  for 
lack  of  a  better  term,  we  have  supposed  to  be  governed  by  the  condition 
which  we  suppose  to  exist,  that  the  velocity  of  efflux  of  the  gases 
from  the  muzzle  becomes  and  remains  sonic  after  the  exit  of  the  pro¬ 
jectile.  There  are  two  cases  to  be  considered.  If  the  velocity  of  the 
projectile  as  It  emerges  is  less  than  the  velocity  of  sound  in  the  gas 
Immediately  behind  It,  we  suppose  that  the  velocity  at  the  muzzle  In¬ 
creases  Instantaneously  to  sonic.  If,  however,  the  velocity  of  the 
pr’ojectlle  Is  greater  than  the  velocity  of  sound  in  the  gas  following 
it,  we  suppose  that  within  the  tube  the  same  Pidduck-IIent  flow  con¬ 
tinues  until  the  velocity  of  that  flow  fails  to  sonic  at  the  muzzle  and 
that  thereafter  the  velocity  of  the  gas  at  the  muzzle  remains  sonic. 
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In  either  ease  once  the  velocity  at  the  muzzle  iB  sonic  the  outflow 
from  the  muzzle  is  more  rapid  than  the  formulas  describing  the 
Pidduck-Kent  flow  would  indicate.  The  resulting  disturbance  of  the 
Pidduek-Kent  flow  is  propagated  through  that  flow  at  sonic  velocity. 
The  equation  of  the  head  of  this  expansion  wave  (or  of  any  character¬ 
istic  in  the  Pidduck-Kent  flow)  may  be  expressed  as  an  elementary  re¬ 
lation  between  quantities  found  by  quadrature  in  the  preceding 
discussion  of  the  Pidduck-Kent  flow.  The  breech  pressure  is  not 
affected  by  this  expansion  wave  until  the  head  of  this  wave  reaches 
the  breech.  In  the  meantime,  the  breech  pressure  is  to  be  computed 
from  the  Pidduck-Kent  formulas. 

To  describe  approximately  the  flow  in  the  expansion  wave  we  have 
used  a  method  which  we  judge  to  be  a  distinct  improvement  over  earlier 
treatments.  It  is  a  flow  in  which  the  Mach  number  remains  constant  at 
each  position,  the  velocity,  pressure  and  density  being  each  a  product 
of  a  function  of  position  by  a  function  of  time.  In  particular  the 
pressure  at  any  point  has  the  same  time  dependence  as  has  been  found 
in  the  well  known  treatment  of  the  emptying  of  a  large  reservoir  of  gas 
through  a  small  orifice  by  assuming  quasi- steady  flow  in  the  orifice 
(see,  for  example,  Aokeret  [l]  ).  The  formulas  so  obtained  have  been 
applied  successfully  to  the  emptying  of  gun  tubes  (see  the  historical 
account  in  Corner  [20  ,  also  Vinti  [ll)  and  Lockett  [73  )  but  not 
happily,  since,  lacking  the  partial  justification  we  give,  here,  the 
effect  of  the  motion  of  the  gas  in  the  tube  could  not  be  assessed. 

The  flow  is  supposed  one-dimensional  and  homentropic.  We  suppose 
the  pressure  p  and  density  p  to  be  connected  by  the  relation  (1.3). 

The  effects  of  friction  and  of  the  motion  of  the  tube  in  recoil  are 
neglected.  Accepting  the  limitations  imposed  by  these  assumptions, 
the  shortcomings  of  the  present  treatment  lie  (l)  in  the  possible 
failure  of  the  Pidduck-Kent  flow  to  provide  a  sufficiently  good 
approximation  to  the  flow  behind  the  projectile,  and  (2)  probably 
chiefly  in  that  the  flow  proposed  as  a  description  of  the  expansion 
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wave  does  not  satisfy  the  conditions  imposed  hy  the  conservation  re¬ 
lations  at  the  head  of  the  wave.  There  is  available  one  parameter 
whose  value  is  to  "be  chosen  to  fit  the  expansion  wave  to  the  Pidduck- 
Kent  flow,  and  while  we  surmise  that  a  proper  choice  would  yield  a  good 
asymptotic  description  of  the  flow  (i.e.  after  a  long  time,  hut  before 
the  time  the  flow  in  the  tube  is  affected  by  atmospheric  pressure  out¬ 
side)  we  have  not  found  a  truly  satisfactory  way  to  determine  the  value 
of  this  parameter.  For  the  purpose  of  estimating  the  behavior  of  the 
breech  pressure  we  have  determined  the  value  of  the  parameter  by  im¬ 
posing  the  condition  that  the  "breech  pressure  be  continuous.  A  nu¬ 
merical  example  is  included  based  on  experiments  performed,  by  Lockett 
[7]  In  these  experiments,  Lockett  measured  the  variation  of  "breech 
pressure  with  time  from  the  instant  of  shot  ejection  until  about  50 
milliseconds  afterwards.  He  then  compared  these  experimental  data  with 
values  predicted  by  two  methods.  In  this  report  we  reproduce  the  re¬ 
sults  of  Lockett,  both  experimental  and  theoretical,  and  adjoin  the 
results  of  our  computations.  These  results  are  shown  in  figures  3  and 
4.  These  are  the  only  pressure  records  known  to  us  for  which  careful 
attention  has  been  paid  to  the  accuracy  of  the  low  pressures  after  shot 
ejection  and  which  also  satisfy  consistency  requirements.  Honesty  re¬ 
quires  the  comment  that  although  our  results  look  very  good  indeed  in 
these  cases,  the  method  cannot  be  expected  to  be  as  good  in  other 
cases,  in  particular  for  low  velocity  guns. 

Figure  5  shows  the  predicted  behavior  of  the  pressure  in  the  tube 
in  one  case.  The  inconsistency  appears  there  as  the  discontinuity  at 
the  head  of  the  expansion  wave.  It  is  clear  that  the  approximation 
is  not  good  in  the  early  stages  of  the  expansion  wave.  For  problems 
In  which  this  portion  of  the  flow  is  important  ve  have  employed  graphi¬ 
cal  integration  using  the  characteristics.  We  hope  to  present  and 
discuss  this  elsewhere.  Let  it  suffice  to  anticipate  that  discussion 
only  to  the  extent  of  remarking  that  the  graphical  calculations  show, 
in  the  cases  we  have  tried,  that  when  the  head  of  the  expansion  wave 
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has  been  once  reflected  at  the  breech,  has  traversed  the  tube  again 
and  been  svept  out  the  muzzle,  the  remaining  flow  in  the  tube  fits  the 
’’expansion  wave”  treatment  presented  here  quite  closely. 


1.  The  Pidduck-Kent  Plow 

The  flow  discovered  by  Pidduck  C8]  and  discussed  by  Kent  [3  and 
Vinti  [12,  133  is  best  described  with  Lagrangian  coordinates.  As  co¬ 
ordinate  to  designate  any  section  of  the  fluid  ve  shall  use  the  mass 
i|r  of  gas,  per  unit  area  of  bore,  behind  it.  Thus  i|f  =  0  is  the  coordi¬ 
nate  of  the  section  of  the  gas  which  is  (and  remains)  in  contact  with 
the  breech.  For  the  gas  at  the  base  of  the  projectile  the  value  of  \[r 
is  the  total  mass  c/g  of  the  charge  divided  by  the  area  A  of  the  bore: 


«  C/Ag. 

Let  the  distance  of  the  section  from  the  breech  at  time  t  be 
x  =  x  (i(r,t).  The  density  is  then 


■  ’  ft 


(1.1) 


If  the  pressure  is  p  =  p  (ty,t),  the  acceleration  of  the  gas  is  given 


5  X  _  Sg 

at.2  ‘  ’ 


(1.2) 


It  Is  supposed  further  that  the  pressure  and  density  are  connected  by 
the  relation 


p  =  k  (i  -  n)“7 


(1-3) 


where  the  constants  q  and  7,  the  covolume  and  the  polytropic  exponent, 
as  they  are  called,  of  the  gas  depend  only  on  the  composition  of  the 
propellant,  and  may  be  considered  known.  The  polytropic  constant  k, 
on  the  other  hand,  specifies  the  particular  adiabatic  curve  followed 
during  the  expansion  of  the  propellant  gas  after  "all  burnt" . 
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Combining  equations  (l.l),  (1.2),  (1.3),  ve  obtain  the  differ¬ 
ential  equation; 

xtt  =~k  ^ 

With  the  free  volume  (as  it  is  called) 


Z  =  X  -  TJ  <|f 
=  Z  (y',t) 

as  independent  variable,  equatioli  (1.4)  takes  the  simpler  form 


ztt  =“k 


(1.5) 


To  solve  this  equation  by  the  method  of  separation  of  variables, 
we  introduce  new  independent  variables,  0  and  w,  and  two  functions  $ 
and  T  by  setting 

0  =  ,  (1.6) 


w  =  et  ,  (1.7) 

z  -  a  0(0)  T  (w),  (1.8) 


where  a,  £  and  €  are  constants  still  to  be  determined.  We  adjust  these 
parameters  to  simplify  the  equation  resulting  from  (1.5)*  Substituting 
(1.6),  (1.7),  (1.8)  into  (1.5)  showG  that  the  relation 


1 


-7+1 

P 


(1-9) 


leads  to 

4.'1  --V5"!"  - 

d  7-1 

where  primes  denote  differentiation  with  respect  to  0  and  w,  res¬ 
pectively.  The  common  value  of  both  sides  of  this  equation  must  be 
constant,  and  without  loss  of  generality  ve  may  assume  this  constant 
to  be  unity.  This  yields  two  ordinary  differential  equations  of  the 
second  order. 
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(7-1)  ($')'7  $"  =2$$' 

(7-l)T"7T'  =  2T’  T" 

From  the  definition,  (0)  =  0.  Designate  as  t  =  0  the  instant 
at  which  the  gas  is  at  rest  so  that  T1  (0)  »  0.  We  set 

1 

$>'(0)  =  1  and  T(0)  =  1  1  (1.10) 

as  two  further  conditions  which  serve  to  specify  completely  the 
functions  <I>  and  T.  Other  non-zero Ichoices  lead  merely  to  a  different 
adjustment  of  numerical  factors  in  <$  and  T  and  their  arguments .  Inte¬ 


grating,  we  have  then 

1  -  T~7+1  =  (T1)2  ,  (1.11) 

1-  ($V7  +  1  =<t>2  ,  (1^12) 

and  the  functions  <J>(0)  and  T  (w)  are  given  by 

— i— 

0  =  J  (l  -  |2)  7-1  d  4  ,  (1.13) 

w  =  J* (1  -  i-’,  +  1)-1/2  45  (1.14) 


Charts,  figures  1  and  2,  furnish  means  for  evaluating  these 
functions.  For  smaller  values  of  and  0,  expansion  of  one  as  a  power 
series  in  the  other  is  convenient.  For  much  greater  precision  the 
tables  of  Vinti  and  his  collaborators  [12,  lj]  may  be  used. 

The  parameters  k,  (3,  oc}  and  e  must  be  evaluated  to  suit  the  indi¬ 
vidual  case.  In  the  type  of  problem  under  consideration  the  muzzle 
velocity  and  the  weight  W  of  the  projectile  are  known,  as  are  the 
we i glit  C  and  chemical  energy  of  the  propelling  gas  and  its  thermo¬ 
dynamic  properties,  7  and  n .  The  values  of  the  parameters  can  be  ob¬ 
tained  from  these  as  is  shown  below.  It  must  be  remembered  however 
that  the  solution  has  r.ot  been  fitted  to  whatever  initial  conditions 
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existed  in  the  gas  when  the  solid  propellant  was  all  hurned.  The  ade¬ 
quacy  of  this  treatment  depends  on  the  conjecture  that  the  asymptotic 
behavior  of  the  solution  is  independent  of  the  initial  conditions  (and, 
of  course,  from  another  point  of  view,  it  depends  on  the  experimental 
verification) . 

Let  the  subscript  p  refer  to  the  base  of  the  projectile.  From 
(1.6)  we  have 

0 

0  =  -  gA  0 _/C  (1.15) 

Y  p  p 

where  A  is  the  area  of  the  cross  section  of  the  bore. 


The  force  on  the  base  of  the  projectile  determines  its  accelera' 
tion  in  the  fashion: 

t  » 


"tt 


where  (1.2)  has  been  used  and  W  is  the  weight  of  the  projectile, 
adjusted  to  account  for  friction.  Inserting  (1.12)  and  (1.15)  there 
results  after  some  reduction, 


_ 7_ 

0  $  (1  -  <£2  )  7  “  1 
P  P  P 


7-  I  C 
2  y  W' 


(1.16) 


which  may  be  solved  simultaneously  with  (1.15)  for  0  and  <I>  as 


functions  of 

W 


Then  p  may  be  found  from  (1.15). 


Let  the  subscript  e  refer  to  the  install  the  base  of  the  pro¬ 
jectile  emerges  from  the  tube.  The  quantities  T^,  T  given  in  (l,l8) 

below  are  found  by  counting  up  the  energy.  The  total  energy  released 

FC 

by  the  propellant  is  ^ - —  ,  where  F  is  the  specific  force  of  the 

propellant.  The  kinetic  energy  of  the  projectile  is  1  .  .  y,v2and  the 

2  g 

ratio  of  these  is  the  efficiency. 


9  * 


7-1 


gFC 


(1.17) 
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and  0  may  "be  computed.  The  internal  energy  of  the  gas  is 


P  P 

JL .  f  pdz=AjL  f  ,-r  dz  =  ^ 

'-ij  *  7-W  t 


A 


7 

7=1 


~o  O 

The  kinetic  energy  of  the  gas  is 


P7 


i  A 


d  if 


^IAoA2  t-7+1  /*  (i_<j>2)  d<J>. 

*<9 


op..«£*’8  7-1  44 


2{S 


We  may  write 

(Total  energy)  =  |  Kinetic  energy  of\+[  Kinetic  energy) +  (  Internal  energy 
\  the  projectile  /  \  of  the  gas  /  \  of  the  gas 

If  the  preceding  expressions  are  inserted  and  some  reduction  carried  out 
there  results 


e(i+-£T.) 

e  37-1  w«<£  2  ^ 


(1.18) 


whence  we  find  T  by  means  of  (l.ll).  The  total  volume  of  the  gun  is, 
e 

from  (1.8) 

^  -  +  ct&  T  A  ,  (1.19) 

g  p  e  * 

This  furnishes  the  value  of  a,  Equation  (i.8)  shows  also  that 


V  =ae$  T ' 
m  p  e 


(1.20) 


which  yields  e.  Finally,  (1.9)  gives  the  value  of  k. 


For  convenience  of  reference,  we  collect  hex e  expressions  for  x^.,p 
and  c  as  functions  of  the  Eulerian  coordinates  x  and  t. 

We  first  calculate^  and  T  as  functions  of  x  and  t.  To  this  end 
we  first  find  w  by  (1.7)  and  T  by.  (l,l4).  Then  we  calculate  0  and  (£> 
by  using  (1.8)  and  (1.13)  simultaneously.  Equation  (1.8)  gives  us 

xt  »  a  e  T1 ,  (1.21) 
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while  (1.1)  and  (1.8)  yield 


p  =  [a  +  tj]  “1 

(1,22) 

The  pressure  p  is,  by  (1.3), 

p  «  k  (a  p^T)"7 

(1.23) 

In  particular,  at  the  breech  <3E>  =  1  and 

^breech  “  k  ^  ^ 

(1.24) 

Finally,  the  sonic  velocity  c  is  found  "by  means  of  the  relation 

2  dp 
c  =t 

dp 

to  "be 

—  7-1 

0  -  (7*)'1/2  (apfc)  2  [l  +  g  J1jT  ]  ■  (1.25) 

If  the  muzzle  velocity  is  less  than  sonic  velocity,  an  expansion 
wave  starts  to  travel  from  the  muzzle  toward  the  "breech  at  the  instant 
of  shot  ejection.  However,  if  the  muzzle  velocity  is  greater  than  the 
sonic  velocity,  there  will  "be  a  time  interval  during  which  the  propellant 
gas  velocity  at  the  muzzle  ctropp  to  sonic.  We  shall  be  interested  in 
calculating  the  duration  of' this  interval.  Let<&  and  T  denote  the 
values  of  3>  and  T  respectively  at  the  muzzle  and  at  the  instant  when 
x.  =  c.  From  (l.2l)  and  (1.25)  we  obtain  one  of  the  conditions  to  be 

v 

satisfied  by<[>  and  T  : 

8  S 

.■Zzi 

1>8  *i  (•',*.)  2  [l  ]  •  &■*> 

The  second  condition  to  be  satisfied  by  <I>  and  T  is  based  on  the 

s  s 

fact  that  the  expansion  wave  starts  at  the  muzzle.  Hence  by  (1.8): 

x  =  3,  0  +  T  ,  (1.27) 

m  p  ss 
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vhere  is  the  equivalent  nuzzle  distance.  Solving  these  two  equations 

fa imiu-taneoua ly  gives  us  the  desired  values  of  $  and  T  , 

s  s 

The  head  of  the  expansion  wave,  the  initial  disturbance  of  the 
Pidduck  flow,  is  propagated  at  sound  speed  with  respect  to  the  gas, 
that  is,  along  a  characteristic  o'  that  flow.  See,  for  example, 
Courant-Friedricks  [5]  ,  p.  106.  Its  progress  is  therefore  governed 
by  the  equation 

d\lf 

at  =  -  pc  * 

If  this  relation  is  expressed  in  terms  of  the  variables  <t>  and  T*  it 
can  readily  be  integrated.  The  result  is 

sin”1  =  -sin-1  T’  +  constant.  (1.28) 


The  constant  may  be  evaluated  by  substituting  the  values  $  and  T1  . 

s  s 

In  particular,  at  the  time  at  which  the  expansion  wave  arrives  at  the 
breech 

T*  =  T£  =  sin  (sin"1  T^  +  sln"1^)  .  (1.29) 

This  should  be  an  improvement  on  a  method  used  hy  Kent  [6]  to  estimate 
the  time  of  the  arrival  of  the  expansion  wave  at  the  breech. 

2.  The  Expansion  Wave 

The  differential  equations  (2.3)  and  (2.4)  below,  for  nearly  one¬ 
dimensional  compressible  flow  in  a  tube  of  variable  cross  section, 
have  solutions  in  which  the  velocity  and  the  speed  of  sound  are. 


respectively. 


and 


t- 1 


C 


t  - 1 


0 


(2.1) 


(2.2) 


where  X(x)  and  L(x)  denote  functions  depending  only  upon  the  distance 
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x  along  the  gun  tube,  t  being  a  suitably  chosen  constant.  A  solution 

in  which  u  vanishes  at  the  breech  and  u  =  c  at  the  muzzle  is  used  to 

describe  the  later  stage  of  the  emptying  of  a  gun  tube.  In  the  case 

of  constant  cross-sectional  area  it  is  found  to  be  possible  to  express 

the  functions  X  and  L  parametrically  in  terms  of  quadratures. 

« 

A  similar  treatment  may  be  applied  to  the  emptying  of  a  vessel 
of  any  shape  through  an  orifice  small  enough  that  we  can  accept  an 
arbitrary  specification  of  the  surface  over  which  the  velocity  is  to 
be  sonic.  We  have  not  investigated  the  possibility  of  proving  the 
conjecture  that  the  flow  from  a  tube  into  a  vacuum  has  asymptotically 
the  character  described  here.  It  would  appear  to  be  somewhat  easier  to 
Investigate  than  the  corresponding  conjecture  about  the  asymptotic 
character  of  the  flow  behind  a  projectile.  For  the  emptying  of  a  vessel 
of  a  shape  other  than  a  tube  the  problem  would  include  inferring  the 
shape  of  the  surface  over  which  the  velocity  is  sonic. 

The  differential  equations  which  govern  the  almost  one-dimensional 
flow  of  a  non-viscous  gas  along  a  tube  of  slightly  and  slowly  varying 
cross-section  A  =  A(x)  are,  in  Eulerlan  form, 

A  pt  +  (pAu)x  =  0  ,  (2.5) 

ut  +  “  “x  ‘  ^  '  (2-4) 

where  u  =  u  (x,t)  is  the  velocity  of  the  gas  and  the  other  letters 
represent  the  same  physical  quantities  as  before.  It  will  be  supposed 
that  the  pressure  and  density  are  connected  by  the  relation 

P  =  k  P7  (2.5) 

It  is  convenient  to  introduce  the  sonic  velocity  c  into  equations 
(2.3)  and  (2.4)  instead  of  the  density.  We  have 

c2  =  dp/dp  =  ?\kp7  ~  1  **  7klj^ 


Zli 

P  7  (2.6) 
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y  _  1  y  -  1  A* 

°t  +  u  cx  *  -  Hr-  cu  tt 


u  +  u  u  =  - - -  c  c 

t  x  7-1  x 


(2.8) 


Introducing  (2.1)  and  (2.2)  there  results 
s(L2  -  X2)  X'  *  L2  -  bX2  -  sA 


(L2  -  X2)  L*  =  -  (1-s)  LX  +  SIX2  A- 


(2.9) 


where 


B  - 


7-1 


For  a  given  A(x)  the  system  (2.9)  could  be  integrated  numerically. 
Problems  such  as  we  are  here  concerned  with  would  have  two  point 
boundary  conditions  which  might  cause  a  bit  of  trouble.  We  have  con¬ 
sidered  only  the  case  in  which  A  is  constant.  A  similar  treatment 
could  be  given  for  the  case  in  which  LA' /A  is  constant.  If  this 

constant  has  the  value  -  ,  the  right  hand  sides  of  (2.9)vanish 

s 

at  X  a  L  and  no  singularity  appears. 


When  A'  is  Identically  zero  an  integral  of  (2.9)  can  be  found 
and  the  solution  can  be  expressed  parametrically  in  terms  of  a 
quadrature.  It  is  convenient  to  introduce  as  parameter  the  Mach 
number. 


M  -  M(x)  »  X/L. 
Equations  (2.9)  mey  then  be  written 

2 

v,  _  1  -  B  M  , 

A  3  - -n 

s(l  -  M  ) 

_  -(1  -  s)M 
L  * 


(2.10) 


(2.11) 
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Also 


and 


Integrating, 


LMr  =  X'  -  ML  * 
1  -  f^M2 


s 


(l“  M2) 


L1  1  dL  -  a (l  -  a)  M 
LM'  “  L  dM  “  “ 


,  2„2 
1-  s  M 


1  -  s 


L  =  L  (1  -  a2V?)  2s 


(2.12) 


(2.13) 


•where  L  is  constant.  Then 
o 


1-  8 


=  ML  -  L  M  (l-  a2V?)  2s 


(2.14) 


Inserting  (2.13)  in  (2.12)  and  integrating,  there  results,  since 
M  =  o  at  x  »  o, 


M 

x  *  sLo  J  (1-  M2)  (1-  e^M2) 


1-  3s 
2s 


dM 


(2.15) 


The  value  of  Lq  Is  determined  by  requiring  that  M  =  1  at  the  muzzle, 
x  =  b,  say. 

The  integral  (2.15)  can  he  expressed  in  terms  of  elementary 

functions  when  —  Is  an  integer,  7  »  5/3,  6/4,  7/5,  etc.  With  a  little 
s  g 

manipulation  it  can  he  evaluated  also  from  Chart  I.  Since  s  is  small, 
a  satisfactory  evaluation  is  obtained  hy  expanding  the  second  factor 
of  the  integrand  in  powers  of  M2.  Errors  of  only  a  fraction  of  one 
percent  are  made  hy  omitting  the  second  factor.  Then 

*  -  <M-  y> 
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and  for  M  =  1, 

x  =  b  =  §  s  L  (2.17) 

J  o 

\ 

Expanding  (2.15)  similarly  ve  have,  approximately, 

L  =  Lq  (1-  m2)  (2.18) 

and 

X  =  MLp  (1-  -s.  -fe-g?.  m2)  (2.19) 

In  this  approximation  the  character  of  the  singularity  at  the  muzzle 
is  unchanged.  Consideration  of  the  first  neglected  terms  indicates 
that  the  maximum  error  in  the  pressure  introduced  In  our  examples  by 
the  use  of  (2.I5)  instead  of  (2.16)  would  be  at  most  about  certainly 
negligible  in  the  presence  of  other  crudities  of  the  treatment.  Actu¬ 
ally  the  plotted  values  were  computed  by  a  different  approximation 
subject  to  errors  of  the  same  size. 

3.  Sample  Calculations 

Lockett  [7]  has  measured  pressures  at  the  breech  of  the  3"/70 
after  the  exit  of  the  projectile  and  to  the  time  when  the  pressure  has 
dropped  to  perhaps  50  psi.  He  has  also  computed  estimates  of  the  pres¬ 
sure  by  the  schemes  which  have  been  proposed  by  Rateau. 

We  shall  now  carry  out  computations  of  the  breech  pressures  and 
of  pressure  profiles  by  the  method  proposed  above  and  compare  the  results 
with  the  pressures  measured  and  computed  by  Lockett.  The  results  of 
these  computations  are  shown  in  figures  3  and  1,  together  with  the 
pressures  measured  and  computed  by  Lockett. 

We  shall  now  explain  the  details  of  our  computations  for  rounds 
No.  1U3-I5I.  We  begin  by  listing  the  data  given  by  Lockett,  converting 
them  to  the  inch- pound  (force) -second  system.  The  charge  is  described 

y  -  ratio  of  specific  heats  =  1.3# 
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r\  »  covolume  =  1.0 6  cm^/s  (in  volume / weight  units) 

^  -1  —2 

=  11,500  in  lb”  sec"  (in  volume/mass  units), 

F  e  force  of  the  powder  =52.9  long  tons/sq.  in. 

POY  g/  cm 

=  _3.20_*  10^  inches, 

C  =  weight  of  charge  a  iq  lbs,  2  oz,  8  drams 
=  10.2  lbs. 

The  gun  is  described  by: 

A  =  cross-sectional  area  of  the  bore,  «  7*°7  sq.  in. 

W  a  weight  of  the  projectile,  «  15  lbs. 

Vm=  muzzle  velocity,  =  3370  ft/sec  =  40,400  in. /sec 
(this  is  an  average  of  the  five  measurements  ob¬ 
tained  for  rounds  143,  14-5,  147,  149  and  15-0. 
total  volume  of  the  gun,  consisting  of  the  volume 
of  the  gun  barrel  and  of  the  volume  of  the 
chamber,  =  1718  cu.  in. 

Note:  U.  atf.  +  AX  ,  where  U  ,  is  the  chamber  volume,  A  the  bore 
t  eh  nr  ch  7 

area  and  X^  the  travel  cf  the  projectile. 

We  shall  now  compute  the  constants  occurring  in  the  theory  of 
the  Pidduck-Kent  flow,  i.e.  a,  fJ,  and  e.  These  constants  will  enable 
us  to  compute  the  pressure,  density,  propellant  gas  velocity  and  sonic 
velocity  in  that  part  of  the  flow  which  is  presumed  to  be  governed  by 
the  Pidduak-Kent  theory. 

We  begin  by  calculating  the  values  of  <$>  and  0  at  the  base  of  the 
projectile.  To  this  end  we  solve  (1,13)  and  (l.l6)  simultaneously. 

The  symbol  W'  of  (l.l6)  designates  the  projectile  weight  after  it  has 


CONFIDENTIAL 


CONFIDENTIAL 


"been  adjusted  to  account  for  energy  losses  due  to  Inertial  and  dissi¬ 
pative  resistance.  We  have  found  it  practical  to  multiply  W  by  1.05 
to  obtain  W1 ,  This  yields  W1  =  15.3  lbs.  In  order  to  solve  (1.15) 
and  (1.16)  simultaneously,  we  may  proceed  as  follows:  we  first  guess 


at  a  value  for  <t>  ,  say  <$> 

P  P 


(1) 


jLzll] 

L  2y  W'J 


1  1/2 


We  use  chart  1  to 


find  the  value  of  0  corresponding  to  $  ,  say  .  We  then  cal¬ 
culate  the  left-hand  side' of  (l.l6),  using  0  =  0p  and  <J>  = 


(!) 

By  linear 


We  then  guess  a  second  value  for  <£>  and  repeat  this  process 
interpolation  ve  then  find  a  third  guess  for  ,  and  we  iterate  in 

this  manner  as  often  as  may  be  necessary.  This  yields  p  =  0.246  and 

0p  =  0.250. 

We  next  find  the  constant  3  by  means  of  equation  (1.15) 

P  =  61.9  in5  lb-1  see"2  . 


To  compute  the  constants  a  and  e,  we  begin  by  finding  the  thermo¬ 
dynamic  efficiency  8  at  the  Instant  of  shot  ejection,  by  means  of 
(1.17) s  0G  -  0.501.  Next,  ve  find  T*  at  the  instant  of  shot  ejection 
by  means  of  (l.l8):  T*  =  0.602.  It  follows  from  (l.ll)  that 
Te  =  4.47. 


We  are  now  ready  to  calculate  the  constants  a,  e  and  k.  Equation 
(1.19)  shows  that  a  -  185  inches,  and  equation  (1.20)  yields 


e  =  I50O  sec 
k  =  2.12  •  1010 


Finally  equation  (1.9)  furnishes 


lb 


-7+1 


-27 

sec 


These  constants  enable  us  to  calculate  the  propellant,  gas  velocity 
x^.,  the  density  p,  the  pressure  p  and  the  sonic  velocity  c  as  functions 
of  the  distance  x  and  the  time  t  during  the  Pidduck-Kent  flow  [cf . 
(1.21)-(1.25)]  . 

However,  the  Pidduck-Kent  flow  ultimately  ceases  when  the  rare¬ 
faction  wave,  starting  at  the  muzzle,  has  reached  the  breech.  Hence, 
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in  order  to  determine  the  time- space  interval  during  which  the  Pidduck- 
Kent  flow  obtains,  we  shall  now  investigate  the  progress  of  the  ex-  . 
pansion  wave. 

We  begin  by  calculating  the  time  when  the  expansion  wave  starts 
at  the  muzzle.  If  the  muzzle  velocity  is  sub- sonic,  the  expansion 
wave  will  start  at  the  muzzle  at  the  instant  of  shot  ejection.  How¬ 
ever,  if  the  muzzle  velocity  is  super- sonic,  the  expansion  wave  will 
not  start  upstream  until  the  velocity  at  the  muzzle  of  the  escaping 
powder  gases  has  dropped  to  sonic.  To  decide  which  case  we  are  dealing 
with,  we  shall  calculate  the  sonic  velocity  at  the  muzzle  at  the  instant 
of  shot  ejection.  For  this  purpose  we  use  (I.25),  which  yields  cg  - 
37 >500  in. /sec.  Comparing  this  with  the  given  muzzle  velocity 

=  40,400  in. /sec,  we  see  that  the  muzzle  velocity  in  the  present 
numerical  case  under  consideration  is  supersonic. 

To  find  T  ,  which  is  the  value  of  T  characterizing  the  instant 
s 

when  the  expansion  wave  starts  at  the  muzzle,  we  solve  (1.26)  and  (1.27) 

• 

simultaneously.  To  this  end  we  guess  at  a  value  for  ,  say  <J>  v  1 , 

Q  3 

and  use  chart  1  to  find  the  corresponding  value  of  0,  say  0^ .  Using 

3 

(1.27),  we  find  the  corresponding  value  of  T,  say  T  ^ .  This  enables 

s 

us  to  find  tb“  right-hand  side  of  (1,26).  A  new  estimate  for  <$,  say 
<3>  is  now  obtained  by  dividing  this  right-hand  side  of  (1.26)  by 

3 

T'  This  process  is  repeated  as  often  as  necessary.  We  obtain 

3 

$  =  0.210,  T  =  5.37  and  T'  =  O.629. 

3  7  3  3 

Having  started  at  the  muzzle,  the  expansion  wave  now  travels 
towards  the  breech.  The  formulae  relating  to  the  Pidduck-Kent  flow 
still  remain  applicable  to  that  portion  of  the  propellant  gas  which  is 
included  between  the  head  of  the  expansion  wave  and  the  breech.  How¬ 
ever,  this  applicability  will  cease,  naturally,  when  the  expansion 
wave  reaches  the  breech.  We  therefore  calculate  next  T^,  which  is 
the  value  of  T  characterizing  the  instant  when  the  expansion  wave 
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arrives  at  the  "breech.  To  this  end  we  use  (1.29),  which  can  "be  put  into 
the  following  more  convenient  form,  which  is  recommended  when  trigo¬ 
nometric  tables  are  not  available: 

2  l/2  2 

T*  »  4>  (l  -T*  )  +  T'  (l-  $  ) 

b  s'  s  '  s'  8  • 

This  furnishes  T^  =  0.779-  Hence,  "by  (1.11),  *  22.4. 

We  have  thus  characterized,  "by  means  of  the  corresponding  values 

of  T,  three  instants:  shot  ejection,  start  of  the  expansion  wave  at 

the  muzzle,  and  arrival  of  the  expansion  wave  at  the  breech,  viz.,  we 

have  found  T  =  4.47,  T  =  5.57  and  T,  =  22.4. 
e  s  b 

To  find  the  values  of  the  time  t  at  these  three  instants,  we  use 
chart  2.  Let  us  consider  the  instant  of  shot  ejection.  Here 

p 

T1  =  O.362.  Using  chart  2,  we  find  w(T  )  =  8.64.  Finally  we  vise 

C  6 

equation  (1.7)  to  calculate  t  =  5-7®  msec.  Similarly,  we  obtain 

©  • 

t  *  6.76  msec,  t^  =  22.4  msec,  all  referred  to  the  instant  t  =  0 
corresponding  to  T  =  1,  For  practical  purposes,  however,  it  is  more 
convenient  to  refer  events  to  the  instant  of  shot  ejection.  We 
therefore  form  the  differences,  to  find  t  -  t  =  0.9^5  msec, 
t^-  tg  ”  16,6  msec.  Thus  the  expansion  wave  starts  0.9^5  msec  after 
shot  ejection,  and  arrives  at  the  breech  16.6  msec  after  shot  ejection. 

We  are  now  ready  to  discuss  the  computation  of  pressure  profiles 
whenever  the  formulae  for  the  Pidduck-Kent  flow  are  applicable. 

1 

In  the  present  case  it  was  decided  to  calculate  the  pressure  at 
4  msec  intervals  from  0  to  16  msec  after  shot  ejection;  and  at  60  inch 
intervals  along  the  gun  tube  from  0  to  180  inches  (measured  from  the 
breech),  as  well  as  at  the  muzzle,  which  corresponds  to  U^/A  =  243 
inches  from  the  breech. 

We  begin  by  calculating  the  progress  of  the  expansion  wave  along 
the  gun  barrel  from  the  muzzle  to  the  breech.  To  this  end  we  first 
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calculate  2  corresponding  to  each  instant  of  time  t  for  which  a  pressure 
profile  is  desired.  Let  us  carry  out  a  sample  calculation  for  the 
instant  of  3  msec  after  shot  ejection.  This  instant  is  suitable,  since 

tQ-  tg  <  8  •  10  J  <  t^  -  tg.  Denoting  this  instant  by  tg,  we  have 

tg  *  tg+8  53  ^5*®  ®sec,  and  hence  by  (1.7)  end  (1.14),  Tg  »  12.6.  We 
rewrite  (1.28)  in  the  form 

$=t£{i-t’2)  -t’(i-t^2) 

Then  <X>  g  =  .07^8.  This  is  a  measure  of  the  mass  of  powder  gas  (per 
unit  area  of  the  bore)  included  between  the  breech  and  the  head  of 
the  expansion  wave  at  the  instant  tg.  The  distance  x  from  the  breech 
to  the  head  of  the  expansion  wave  at  this  instant  can  now  be  calculated 
by  equation  (1.8)  yielding  x(tg)  =  186  inches.  We  compute  £{t)  In  this 
manner  for  the  prescribed  values  of  t  satisfying  t  <  t  <  t^. 

We  are  now  ready  to  compute  the  pressure  profiles  for  that  portion 
of  the  flow  which  is  described  by  the  Fidduck-Kent  formula  (1.23).  For 
values  of  t  satisfying  t  <  t  <  t  ,  we  can  find  the  pressure  p  as  a 
function  of  the  distance  x  (measured  from  the  breech)  for  all  values 
of  x  along  the  gun  tube,  i.e.  for  all  x  satisfying  0  <  x  <  IJ^/A.  For 

all  values  of  t  satisfying  t  <  t  <  t,  ,  we  can  find  p(x)  for  all  x 

A  S  —  D 

satisfying  0  <  x  <  x(t) .  But  the  flow  in  the  remainder  of  the  semi¬ 
infinite  strip  cf  the  x  - 1  plane  0  <  x  <  U  /A,  t  <  t  is  not  described 
by  the  Pidduck-Kent  formulae. 

Let  (x,t)  denote  a  point  of  this  semi-infinite  strip  for  which 
the  Pidduck-Kent  formulae  are  applicable.  To  calculate  the  pressure 
p(x,t),  we  begin  by  finding  corresponding  to  the  given  point  (x,t). 

As  an  example,  let  us  consider  the  point  120  inches  from  the  breech  a 
and  8  msec  after  shot  ejection.  Since  120  <  l86  =  ^(tg) ,  the  Pidduck- 
Kent  formulae  are  applicable  at  this  point.  To  find  corresponding 
to  ve  use  (1.8),  whence  $  (x^qj^  =  The  pressure 
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is  now  found  by  equation  (1.25):  -  4200  lbs/sq.  in.  The 

pressure  at  the  instant  tg  is  thus  computed  for  all  desired  values  of 
x,  including  e.g.  the  pressure  at  the  head  of  the  expansion  wave,  by 
setting  x  =  x  (tg).  In  our  example  we  obtain  p  =  4l4Q  psi  for  this 
pressure. 

Pressure  profiles  can  thus  be  calculated  for  all  desired  values 
of  t,  for  the  Pidduclc-Kent  portion  of  our  flow. 

In  particular,  we  need  to  know  the  pressure  at  the  breech  at  the 
instant  the  expansion  wave  arrives  there,  i.e.  at  t^.  Since  in  our 
example  we  have  =  22.4,  (1.24)  shows  that  the  breech  pressure  at 
t^  is  p(o,t^)  =  2010  psi. 

The  formulae  for  the  flow  of  the  second  type  are  assumed  to  apply 
to  that  portion  of  the  flow  which  is  not  described  by  the  formulae  of 
the  Pi&duck-Kent  flow,  i.e.  for  tp<  t<_  t^,  x(t)  <  x<  U^/A,  and  for 
^  <  t,  o  £  x  ^  yA. 

Various  conditions  could  be  imposed  in  order  to  effect  the  deter¬ 
mination  of  t  .  Since  we  are  interested  in  variations  of  the  breech 
o 

pressure  with  time,  we  have  imposed  the  requirement  that  the  breech 
pressure  be  continuous  at  the  instant  of  arrival  of  the  expansion  wave 
at  the  breech,  t  =  t^  =  .0224  sec.  We  had  found  above  that  p(0,t^)  = 
2010  psi  and  the  corresponding  velocity  of  sound  1b  c^  =25,200  in. /sec. 
Since  we  are  neglecting  covolume  excess  from  this  point  on,  the  manner 
in  which  c  should  be  computed  is  indeterminate  to  the  extent  of  about 
2$.  Fortunately  the  effects  of  this,  ultimately  17$,  become  large 
only  at  later  times  when  the  percentage  accuracy  of  our  determination 
and  of  the  measurements  is  equally  reduced.  From  (2,17)  we  find 

5U 

L  =  —  =  =  2430  in. 

o  2s  . 5A 

and  (2.2)  with  x  =  M  =  0  gives 


L 

0 
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from  which  t  =  -.0741  sec. 
o 

Having  determined  tQ,  we  are  ready  to  calculate  the  desired  pressure 
profiles.  Let  us  calculate  e.g.  the  pressure  at  t^Q,  which  is  the 
instant  20  msec  after  shot  ejection,  and  at  xQq  =  90  inches  from  4 
"breech.  We  note  first  that  t  =5.78  msec  implies  t„  =25.8  msec. 

From  (2.16) 


M  - 


.24-7 


and  we  may  solve  for  M  =  .252.  The  pressure  is,  using  (2.2)  and  (2.13) 

_2y 

*V  v  c. 


p  *  p(  QiO  ( 7  )  7  -1 


2y 

-  1(i-s2m2) 


7(3-  y) 
(7- 1)2 

(3.i) 


=  1450  psi, 


when  the  values  above  are  inserted. 

As  a  second  example,  let  us  compute  the  pressure  at  the  head  of 
the  expansion  wave  8  msec  after  shot  ejection.  By  our  previous  cal¬ 
culations  we  had  found  that  x(tg)  =  186  inches.  Hence  by  equations 
(2.16)  and  (3-1)  p[x(tQ),  tgj  =  3770  psi,  as  compared  to  4-140  psi 
calculated  at  the  same  instant  and  at  the  same  point  by  the  Pidduck- 
Kent  formulae. 


Pressure  profiles  can  thus  he  calculated  for  all  values  of  t  >  te, 

either  hy  the  Pidduck-Kent  formulae  or  by  the  formulae  for  the  flow 

of  the  second  type.  The  former  hold  for  t  <  t  <  t  and  all  x,  and 

A  e  s 

for  t  <  t  <  t,  and  0  <  x  <  x(t);  while  the  latter  are  to  be  used  for 
S  D  . 

t^  <  t  <  t^  and  x(t)  <  x  <  Uj./A,  and  for  t  >  and  all  x. 


WILLIAM  C.  TAYLOR 


P v  M  •  /  —  1^)  . . 

PETER  M.  TREUEHFEL3 
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